Astronomy Sz Astrophysics manuscript no. HeatfluxAA 


February 2, 2008 


(DOI: will be inserted by hand later) 





Thermal structure of gas in relaxed clusters of galaxies 

E. Fischer 

Auf der Hoehe 82, 52223 Stolberg, Germany 



O 
O 

(N 
G 



> 

O 
O 

o 



Oh' 
i 1 

O 1 



X 



Received 2.1.2004 / Accepted 

Abstract. Gas clouds under the influence of gravitation in thermodynamic equilibrium cannot be isothermal due to 
the Dufour effect, the energy flux induced by density gradients. In galaxy clusters this effect may be responsible for 
most of the "cooling flows" instead of radiative cooling of the gas. Observations from XMM-Newton and Chandra 
satellites agree well with model calculations, assuming equilibrium between energy flux induced by temperature 
and density gradients, but without radiation loss. 
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1. Introduction 

It has been known for long years that in most clusters of 
galaxies the intracluster plasma is far from being isother- 
mal, exhibiting strong decrease of the plasma temperature 
in the vicinity of individual large galaxies. This effect is 
commonly attributed to "cooling flows", that means, in- 
creasing radiative energy loss with increasing density of 
the gas, which is gravitationally attracted by the galaxies. 
But recent high resolution data from the XMM-Newton 
and Chandra satellites have shown that this cooling flow 
model does not fit to the observations. This has newly 
kindled a vivid discussion on the the origin of the discrep- 
ancy. 

Already in former years the fact was addressed that 
thermal conduction would wash out the temperature gra- 
dients, if the thermal conductivity of the plasma were 
of the magnitude given in the classical book by Spitzer 
(1962). As a possible explanation it has been proposed 
that magnetic fields may reduce the conduction. Several 
authors claim very different reduction factors of the ther- 
mal conductivity (see e.g.Voigt et al. 2002, Markevitch et 
al. 12)03) to fit the data. The existence of magnetic fields 
in the order of 0.1 — 1/xG has been confirmed by various 
observations (see e.g. Dolag et al. 2001). But to produce 
the required reduction, ordered magnetic fields of a still 
higher magnitude would be necessary, as the heat flux re- 
duction works only perpendicular to the magnetic field 
lines and there is not yet a generally agreed model to de- 
cide, to which extent magnetic fields can influence thermal 
conductivity. Narayan and Medvedev 120011 have shown 
that fields fluctuating on a large range of length scales 
can produce only minor changes in the thermal conduc- 



tivity, while others (see Chandran and Maron 2004] and 
references therein) give reduction factors as high as 10. 

One problem with all these reduction factors is, that 
they would affect the temperature profiles preferably at 
high temperatures in the outer boundary of the " cooling 
flow" region, while radiative cooling is most effective near 
the centre. But observations exhibit just the opposite ef- 
fect (Kaastra et al. EffiMJl 

In all these discussions it is implicitly assumed that 
in thermodynamic equilibrium conduction leads to a uni- 
form temperature, though it is well known that density 
gradients may also cause a flux of energy and thus induce 
temperature gradients, an effect known from the text- 
books as "Dufour effect" (see e.g. Hirschfelder et al. 1964). 
Under terrestrial conditions this effect is regarded mostly 
as a second order correction to "normal" conduction, that 
means, to energy flux due to temperature gradients. We 
are more familiar with the "reciprocal" effect, the diffusion 
of particles caused by temperature gradients, known as 
thermal diffusion. But under the condition of cosmic plas- 
mas, where stationary density gradients are maintained 
by gravitational fields, the Dufour effect can considerably 
alter the equilibrium conditions. Even in relaxed systems 
there will be stationary temperature gradients to balance 
the energy flux caused by density gradients. 

Normally in a gas we associate " equilibrium" with con- 
stant pressure, density and temperature. But in a system, 
which is influenced by volume forces such as gravitation, 
the equilibrium state may well exhibit gradients of the 
state variables. Mass flow equilibrium in these systems is 
obtained by the balance between the gravitational force 
and the counteracting pressure gradient, which is related 
to the density gradient by some equation of state. 

But if there is a density gradient, energy equilibrium 
can be obtained only, if there is also a temperature gradi- 
ent, balancing the Dufour effect, to obtain zero net energy 
flux. Below we will derive the equilibrium conditions for a 
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system, which is fully relaxed with respect to flux of mass, 
momentum and energy, within the scope of a simplified ki- 
netic model. 

In the next section we will discuss the equilibrium con- 
ditions in cluster plasmas. In the third section a simpli- 
fied kinetic model is developed to determine the essential 
consequences of the Dufour effect. In section 4 the results 
are compared with observational data from XMM-Newton 
and Chandra, and in section 5 conclusions are drawn with 
respect to the " cooling flow" region as well as for the equi- 
librium conditions of complete galaxy clusters. 

2. Equilibrium relations 

Speaking of equilibrium means that the net fluxes of all 
properties determining the state of a system are locally 
balanced. Complete local equilibrium will of course never 
be established in a system with radiation losses or with ir- 
reversible chemical reactions. But as long as these are not 
the dominant exchange processes, quasi-equilibrium con- 
ditions can be formulated, considering only the dominant 
transport processes. 

To identify the dominant transport processes, it is con- 
venient to compare the time constants for the local change 
of temperature, density or composition due to some loss 
and production process or flux. Common to all discussions 
on the energy balance of the " cooling flow" region of intra- 
cluster gas is, that radiation is regarded as the only energy 
loss process. This loss has to be balanced by energy gain 
from inflowing hot matter, by uptake of gravitational en- 
ergy from inflowing gas, by heat conduction from outer hot 
regions and energy production by supernovae in galaxies 
or active galactic nuclei. 

The Dufour effect adds a further transport process to 
the balance, which acts like an additional cooling, op- 
posite to the heating processes mentioned above. It can 
be regarded as a thermal sedimentation of low energy 
particles into the vicinity of a central gravitating mass. 
Microscopically it is closely connected to the "normal" 
thermal conduction, as both result from the same process: 
the transport of kinetic energy of the plasma particles in 
an inhomogeneous medium. If there is thermal conduc- 
tion, there is always one fraction caused by temperature 
gradients and one caused by density gradients. The rela- 
tion between the two parts will be discussed more in detail 
in the next section. 

To find out, which of the heating or cooling pro- 
cesses preferentially contributes to the balance we com- 
pare the time constants of the corresponding processes. 
In the "classical" cooling flow model only two processes 
are considered: radiation losses and inward flow of hot 
matter. The time constant for radiative cooling is r r = 
(3nkT)/(A(T) n 2 ), where A(T) is the cooling function (see 
e.g.Voigt et al. 12002). In the temperature and density 
range observed in the centre of galaxy clusters typical val- 
ues are kT — 3 keV and n — 10~ 2 cm~ 3 , resulting in a cool- 
ing time r r = 1.2 x 10 17 sec. As this is much longer than 
the collapse time t c = 1 / y 1 pG = 3x 10 16 sec, equilibrium 



can always be reached between gas inflow and radiation 
loss. 

There are, however, two mayor problems with this 
model. It does not fit to the observations, as the increas- 
ing cooling rate should increase with increasing density, 
so that the temperature profile should exhibit a sharp 
dip in the centre, which is not observed. Secondly it ne- 
glects thermal conduction, which has time constants much 
shorter than radiation, so that it could wash out all gra- 
dients established by radiative cooling. If thermal con- 
duction were at the theoretical value given by Spitzer 
UJJjlA = 4.6 x 10~ 7 T 5 / 2 erg cm-^-^ec" 1 , for a cooling 
flow region of r c = 50 kpc the conduction time scale is 
T\ = nkr 2 /X = 10 16 sec, much less than the cooling time 
scale. Even if thermal conduction is reduced by the pres- 
ence of magnetic fields, as proposed by several authors 
(see Narayan and Medvedev 120011 Chandran and Maron 
12777741 . it remains a dominant transport process. 

But as has been mentioned above, this does not im- 
ply that the temperature should be constant. Due to the 
opposite effects of temperature and density gradients also 
under the condition of local thermal equilibrium the tem- 
perature will vary in radial direction. In the next section 
we will derive this equilibrium condition under the as- 
sumption that conduction is dominant and radiation losses 
can be neglected. 

3. Kinetic model 

To determine the equilibrium conditions in a plasma un- 
der the influence of an outer force field such as gravitation, 
we determine the net flux of mass and energy through a 
control area AF under the condition that the particles in 
the gas have a Maxwellian distribution everywhere, but 
with number density and mean energy varying perpendic- 
ular to the control area. The net flux is calculated then 
under the assumption, that any particle crossing the con- 
trol area retains its velocity and direction, which it has 
obtained in the last collision, one mean free path from the 
surface, and in addition is subjected to some acceleration 
from the force field. 

In the energy balance the gain and loss of kinetic en- 
ergy due to the changing gravitational potential will be 
omitted, as it cancels out from the balance exactly, when 
the net particle flux through the control area is zero. As 
we are interested only in the equilibrium condition, the 
absolute magnitude of the energy or particle flux is not 
important, as long as it dominates over other transport, 
production or loss processes, which appears to be justified 
for most of the "cooling flow" region in clusters. 

Denoting the direction perpendicular to the surface as 
z and the angle between this and the direction of particle 
motion as 9, the flux of some property G through the 
surface is 

j = J AF cos d(j+-j-)<m, (1) 

where j + and j~ are the normal components of the flux 
of G carried by particles moving towards the surface 
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in positive and negative z direction under the angle 9. 
The integration has to be taken over half the solid angle 
7t/2 > 9 > 0. Denoting the amount of property G car- 
ried by particles moving with velocity v by g(v), the flux 
components of g(v) are 



= (a(v) - ^r^Acos0)(wcos0 + &-) 
dz v 



(2) 



Introducing eqs.© 
mass and energy flux 

dn 



to into eq.(Q}, we get for the 



ju = 2B a [ xB— + x(2x l - 3)n 



dp 
dz 



2bn 1 
~B~x 



h(x)dx (8) 



JE = 3 a+2 [ x 3 B^+x 3 (2x 2 - 3)nf—^x) h{x)dx (9) 
dz dz p J 



J 



dg(v) 

dz 



A cos 9) (v cos ( 



v 



2+a, 



(3) 



A is the mean free path, which may depend also on the 
velocity. The last term denotes the acceleration of the par- 
ticles during their flight from the last collision to the con- 
trol surface. Integrating eq.Q over the solid angle yields 
the equation 



j = -ttAF Aw 



dg(v) 
dz 



2b\ 



-g(v) 



(4) 



To determine mass and energy transport by particles 
moving with velocity v, g has to be set to 



g M = n(z)mf(v,z)dv, g E = n(z)m— } \v , z)dv (5) 

In case of a neutral gas n(z) and m are the numerical den- 
sity and mass of atoms, f(v, z) is the distribution function 



4 v 2 _« 2 
f(v)dv = —j=-^e ~^dv 
v 71 " P 



(6) 



B(z) being the abbreviation B{z) = y / 2kT(z)/m. To cal- 
culate the quantity dg/dz we also need the derivative 



dz 



d(3 
dz 



2v 4 



/? 4 



Pdz\P 2 " 



f(v)(7) 



While in a neutral gas the parameter b in eq. J3J is the 
normal gravitational acceleration, in a plasma the action 
of the force is somewhat indirect. Energy transport is me- 
diated preferably by electrons. Thus in eqs.© n(z) and m 
are the electron density and mass. But in this case acceler- 
ation by the external force is acting also on the ions, which 
transfer the force to the electrons by Coulomb interaction. 
Thus the parameter b in eqs.© and ® is not the gravi- 
tational acceleration of free moving electrons. But for the 
equilibrium condition of zero transport the absolute mag- 
nitude of b is not relevant, when we only want to know, 
under which conditions the net flux of electrons and the 
flux of energy carried by them are in balance. 

An additional difference between neutral gas and 
plasma results from the fact that the collision cross sec- 
tion between neutral atoms is nearly independent of the 
collision energy, so that the mean free path A does not 
depend on v. In plasmas, where Coulomb interaction is 
dominant, the mean free path increases with v 2 . Thus in 
the further calculations we set A = Xo(v/vo) a with a = 
for neutral gas and a — 2 for a fully ionised plasma. 



with x = v/(3 and h(x) = 2y / 7rrnA_F(Ao/wo ) x 
Integrating these equations over all x, we obtain the con- 
ditions for zero mass and energy flux: 



dn _ dB 
B—F 3+a +n-^[2F 5+a 
dz dz 



3F : 



2bn 



3+aJ 



'-Fu 



(10) 



dz 



2bn 



— F3 +a - (11) 



The abbreviation Ff. stands for the integral F^ = 
x k e~ x dx. The different arguments of F^ in eas. (|l(Jfl 
and (|11|) can be eliminated by the recurrence formula 
Fk+2 = (k + l)/2 x Fk, so that we finally obtain the equi- 
librium conditions for mass and energy transport: 



+ a n dn (2 + a)(l + a) d/3 2bn 

— (3 Tz + 2 n Tz~ — = 



A + adn (4 + a)(3 + a) dB 2bn 
8 h -n— 

2 ' dz 2 dz B 







(12) 



(13) 



It is immediately evident that this is a set of linear 
equations for the gradients dB / dz and dn/ dz. If there are 
no external forces (b — 0), there exists only the trivial 
solution dB / dz = and dn/dz = 0. 

In the presence of volume forces such as gravitation 
eas. H12|l and (|13|l are inhomogeneous and allow a non- 
trivial solution, which constitutes a fixed relation between 
temperature gradient and particle density gradient or, be- 
cause of quasi-neutrality, also between temperature and 
mass density gradient. Eliminating b leads to the relation 



„dn , . dB 

B— + (2a + 5)n--=0 
dz dz 



(14) 



with the solution nB^ 2a+b ^ = const. Replacing the number 
density by the mass density p and B by ^/2kT/m we finally 
find the condition pT 5 / 2 = const, for neutral gas and 
pT 9 / 2 — const, for a fully ionised plasma. 

This result is independent from the actual magnitude 
of the coefficient of thermal conduction, if it is influenced 
by magnetic fields or not. The only conditions, which have 
to be satisfied, are the existence of local thermal equilib- 
rium and the condition, that transport of kinetic energy 
of particles is the dominant transfer mechanism. 



4. Observations 

During the last years several spatially resolved x-ray emis- 
sion measurements of galaxy clusters have become avail- 
able, especially from the Chandra and XMM-Newton 
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satellites. Though temperature determinations are still 
rather poor due to projection effects and necessary back- 
ground corrections, intensity profiles and from these den- 
sity profiles have been obtained with high accuracy. But 
in spite of the remaining error margin the classical cooling 
flow model can definitely be ruled out as an explanation 
of the observed data (see Kaastra et al. l2004f) . Neither the 
radial temperature profiles fit to those expected from this 
model nor do the spectra show the expected intensity of 
metal lines of lower ionisation levels. 
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Fig. 1. kT*p 2 / 9 for consecutive annular channels starting 
from cluster center for the 7 most luminous clusters given 
by Kaastra et al. \W01\ 



Good agreement with observations can, however be ob- 
tained with a model based on thermal conduction with in- 
clusion of the Dufour effect. We have applied the model to 
the data given by Kaastra et al. 2004. In fig^the quantity 
T x p 2 / 9 is plotted for the seven most luminous clusters, for 
which data are given by Kaastra et al. Only the clusters 
with the highest X-ray flux (> 5 x 10~ n erg cm _2 sec _1 , 
0.1-2.5 keV) have been used, as for these clusters the tem- 
perature profiles appear sufficiently accurate. The data 
points in figH refer to temperatures in consecutive an- 
nular channels starting from centre to the boundary of 
the cooling region. (For details see Kaastra 20.04). For all 
seven clusters the constancy of T x p 2 ! 9 over the cooling 
region is evident. This can be regarded as a confirmation 
that thermal conduction with inclusion of the Dufour ef- 
fect is the dominant energy transfer mechanism and that 
radiative cooling is less important. 

Data from Perseus cluster have also been obtained by 
Chandra (see Sanders et al. {2004). Fig 13 shows a com- 
parison of the temperature profiles obtained from direct 
measurements by Chandra and XMM-Newton compared 
with data derived from the corresponding density profiles, 
applying pT 9 / 2 = const. The curves are fitted at the clus- 
ter centre. As can be seen, the curves derived from den- 
sity data agree very well with each other and with the 



direct measurements from XMM-Newton, while the tem- 
peratures from Chandra deviate considerably in the outer 
region. It should be mentioned, however, that the Chandra 
data are obtained by averaging angular resolved profiles 
which exhibit a large spread (> 1.5 keV) in the outer re- 
gion. 
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Fig. 2. Temperature profile of the Perseus cluster. 
Comparison of measured temperature with profiles derived 
from measured density by pT 9 / 2 = const., fitted at the clus- 
ter centre. 

Whether radiative loss can completely neglected or has 
at least to be added as a correction to the temperature pro- 
file resulting from pure heat conduction, cannot be decided 
from the present data. For the model this would require a 
better knowledge of the absolute magnitude of the ther- 
mal conductivity. This requires an accurate theory of the 
reduction effects by magnetic fields and more knowledge 
on the distribution of these fields over the cluster volume. 



5. Conclusions 

We have shown that thermal conduction appears to be the 
dominant energy transport mechanism and that a pure 
conduction model can describe the observed temperature 
profiles, when energy flux induced by density gradients is 
added to the balance. The relative temperature gradients 
are small compared to the associated density gradients. 
But temperature changes by a factor of two or three, as 
they have been found in the central regions of the clusters 
mentioned above and which are reported also for other 
"cooling flow" clusters by several authors (see e.g. Voigt et 
al. I2002[) . can well be understood from the Dufour effect, 
as the gas density may change by a few orders of mag- 
nitude in the vicinity of large galaxies. It appears likely 
that most of the observed temperature drop in the " cool- 
ing flows" around these galaxies is not caused by radiative 
cooling but by the "conductive cooling" associated with 
the density gradients. 
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It should be noted that the Dufour effect is not re- 
stricted to the cores of galaxies, but it should be impor- 
tant also in the hot corona of individual stars, where it can 
help to explain the observed temperature gradient from 
the surface into the interstellar medium. 

Applying the model to complete galaxy clusters leads 
to the result that the plasma temperature should con- 
tinuously rise with increasing distance from the centre. 
Observational data are contradictory in the outer region. 
Some authors claim a temperature drop at the outskirts 
of the observable region, others find constant or increas- 
ing temperature (see Mushotzky 2004 and references given 
there) . Most of the recent papers find increasing temper- 
atures at least out to the point, where the increase of the 
error bars exceedes the expected temperature gradients. 

Also the conduction model should be taken with care 
at the low densities in the outer regions, as the conditions 
of local thermal equilibrium are no longer satisfied due 
to the extremely large mean free path of electrons. But 
it can be stated that, if thermal conductivity were close 
to the Spitzer value, clusters would have cooled off in the 
time since formation, unless the clusters are embedded in 
a very hot but thin intergalactic plasma. The existence of 
a diffuse high energy x-ray background would be a strong 
hint to the existence of such an intergalactic plasma. 

Measurements of the background in the energy range 
above that of the intracluster gas are rare. Observations 
from the HEAO-1 satellite (Boldt ITW|l have found 
a continuous background, which can be interpreted as 
bremsstrahlung in the temperature range of 40 to 50 keV, 
but the spatial resolution of these early data is so poor 
that it cannot be decided, if it is really radiation from 
a continuously distributed plasma or from the superposi- 
tion of individual point sources. As a plausible explanation 
for the formation of a hot continuous background plasma 
one could think of the relativistic matter jets expelled by 
quasars. During the time period of several Gyr this matter 
should have become homogeneously distributed in space. 

It has been shown that the conduction model presented 
in this paper can explain some observations in the central 
regions of the intracluster plasma, but it deserves more 
confirmation by detailed measurements, until the astro- 
physics community will be convinced. Also the extension 
to the outer regions of galaxy clusters would bring addi- 
tional confirmation, but this requires the extension of spa- 
tially resolved measurements into the high energy range, 
as they are planned with the launch of future satellites 
like Astro-E2, which are scheduled for the next years. 
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